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1 Introduction 



The present day standard model of cosmology with a cosmolological constant and pressure- 
less dark matter fits oberservational data. However, despite quite some effort, no particle 
candidate for dark matter has been discovered. On the other hand, the standard model 
has only few challenging models with a finite number of parameters. We propose to use 
an Einstein-Cartan's theory as an alternative. The cosmological model we will discuss has 
the same number of parameters as the standard model. The dark matter density today 
of the standard model is replaced by a mean spin density of baryonic matter today. 

Einstein-Cartan's theory [1] has several popular motivations: 

• Space-time torsion, that Einstein puts to zero from the start, is allowed without 
however modifying the Einstein-Hilbert action. Consequently energy-momentum 
is still the source of space-time curvature with Newton's constant being the cou- 
pling constant and the source of torsion is half-integer spin with the same coupling 
constant. 

• To describe matter with half-integer spin, one must use an orthonormal frame and 
the spin connection u. Consequently the gauge invariance of the Einstein-Hilbert 
action is manifest with the gauge group being the Lorentz group. Furthermore, 
Einstein-Cartan's theory treats this connection as an independent field variable 
besides the metric. Both features look like promising steps towards unification of 
gravity with the standard model of particle physics, which is a gauge theory with a 
connection as independent field variable. 

• Torsion closes the diagram of Figure 1. Let us explain the two arrows labeled 'ge- 
ometry'. One possibility to define and measure curvature is by parallel transporting 
a vector around an infinitesimal geodesic parallelogram. The Riemann tensor eval- 
uated on the two vectors defining the parallelogram is precisely the infinitesimal 
rotation mapping the transported vector after the round trip onto the initial vec- 
tor. Torsion can be defined and measured in a similar fashion. In presence of 
non-vanishing torsion, the infinitesimal geodesic parallelogram does not close. The 
translation from the final to the initial point of the 'parallelogram' defines a tangent 
vector, which is precisely the torsion tensor evaluated on the two vectors defining 
the 'parallelogram'. 

Since a long time physicists have been attracted by the beauty of torsion and there 
is a vast body of literature on the subject including applications in cosmology. Let us 
cite the historical review by Hehl, von der Heyde, Kerlick & Nester [2] (which was the 
second author's first contact with torsion), two recent reviews |3] and a few contributions 
of dynamical torsion in cosmology [H E]- Other results directly related to the present 
work with its non-propagating torsion will be referred to as we go along. 

Let us mention that supergravity, which has been remarkably popular, is an Einstein- 
Cartan theory with matter being a spin | Rarita-Schwinger field. If there torsion is put 
to zero, the theory loses its supersymmetry. 
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Einstein eq. 

energy-momentum > curvature 




translations rotations 




torsion < spin 

Cartan eq. 



Figure 1: Einstein and Cartan's equations, geometric definition of curvature and torsion, 
Noether's theorem 

Let us also mention two recent results concerning torsion in noncommutative geometry: 
(%) Torsion is produced together with curvature when the flat Dirac operator fluctuates 
[6 J in the sense of Connes [7j. 

(ii) Connes' spectral action [7] produces torsion terms together with the Einstein-Hilbert 
action 131. 



2 Einstein-Cartan's theory in a nut shell 

It is convenient to write Einstein-Cartan's theory using an orthonormal frame, un repere 
mobile using Cartan's words. Let x M (// = 0,1,2,3) be a coordinate system on an open 
subset of M 4 and let e a =: e a fJi dx tJ- , a = 0, 1,2,3 be four 1-forms which are orthonormal 
with respect to a given space-time metric g. Following the conventions of reference [9], we 
write a metric connection with respect to this orthonormal frame as a 1-form with values 
in the Lie algebra of the Lorentz group, u a b ='■ uj a b^ dx M . Then (suppressing all wedge 
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symbols) Cartan's two structure equations read: 

R:=du + ±[w,w], (1) 

for the curvature, a 2-form with values in Lie algebra of the Lorentz group, R a b —'■ 

\R a bcd e c e d , and 

T := De = de + ue, (2) 

for the torsion, a vector-valued 2-form, T a =: ^T a bc e b e c . It will be useful to decompose 
the torsion tensor into its three irreducible parts: 

Tabc = A abc + r] ab V c - 7] ac V b + M abc , (3) 

with the completely antisymmetric part A abc := |(T afec + T cab + T bca ), the vector part 
V c ■= \T abc r] ab , and the mixed part M abc characterized by M ahc = -M acb , M abc i] ab = 0, 
and M abc + M cab + M bca = 0. 

We have the two Bianchi identities: 

BR = dR+ [to, R] = 0, (4) 
DT = DDe = Re. (5) 

In these notations, the Einstein-Hilbert action reads 

(R ab +lAe a e b )e c e d e abcd 
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J(R ab ab + 2A)dV, (6) 



where * is the Hodge star of the metric g and 60123 = 1- The energy- momentum current 
is the vector-valued 3-form r a obtained by varying the orthonormal frame in the matter 
Lagrangian: 

C M [e + M- £u[e,uj} =: -f a r a + 0{f). (7) 

The energy-momentum tensor r ab is defined by *r a =: r ab e b . 

Likewise, the spin current is the Lorentz-valued 3-form S ab obtained by varying the 
connection in the matter Lagrangian: 

C M [e,u + X ] -Cu[eM =: -\ X ab S ab + 0( X 2 ). (8) 

The spin tensor S abc is defined by *S ab ='■ S abc e c . 

Now we can derive the Einstein equation by varying the total action with respect to 
the orthonormal frame: 

(R ab +lAe a e b )e d e abcd = -1QtvGt c or equivalently G ab - Arj ab = 87rGr ab , (9) 
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with the Einstein tensor G ab '■= R c a cb — \R CC l cdVab- 

Likewise Cartan's equation is derived by varying the total action with respect to the 
connection: 

T c e d e abcd = -8nGS ab , (10) 

or equivalently: 

A ccl b + 3V a r} bc - 3V b r} ac + M cab = -8nGS abc . (11) 

Unlike curvature, torsion does not propagate: it is non- vanishing only inside matter with 
half-integer spin. 

Note that in presence of general torsion, the Einstein tensor G ab and energy-momentum 
tensor r ab are neither covariantly conserved nor symmetric. Indeed combining Einstein's 
equation with the first Bianchi identity, we have: 

(R ab + Ae a e b ) T d e abcd = -16nG(Dr) c , (12) 

and combining Cartan's equation with the second Bianchi identity, we have: 

Tab - r ba = *(DS) ab . (13) 

In absence of half-integer spin, general relativity is usually written with respect to a 
holonomic frame dx^ and the metric g is given by its metric tensor g^x) with respect to 
the coordinates x M . The holonomic frame will be useful below to compute geodesies and 
invariant connections. The metric tensor reads g^x) = e a ^(x) e b v (x) r\ ab . Traditionally 
the components of the connection with respect to a holonomic frame are denoted by 
T"^ =: T'^dx^, a g£(4) valued 1-form. The link between the components of the same 
connection with respect to the holonomic frame T and with respect to the orthonormal 
frame u is given by the GL{A) gauge transformation with e(x) = e a /1 (a;) G GL(A); 

u = ere" 1 + ede" 1 , (14) 

or with indices: 

(9 

, ,a -pa — 1/3 , ^ -la /i r\ 

w bp — e Q l p^e b + ea ~Q^l e b ' 

The given connection (with or without torsion) is metric if 'angles' and 'lengths' are 
preserved under parallel transport. For the components u of the connection this just 
means that its values u a b are in the Lie algebra of the Lorentz group: oo ab = —uj ba . In 
terms of the holonomic components the metricity of the connection reads: 

9^ ~ I^p,, - V\ x g^ = 0. (16) 



dx x 



In presence of general torsion, there are two different geodesies x lx (p). The first is 
defined by parallel transport of the velocity vector x with the given connection: 

x A + rVi"f = 0, (17) 
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where the over-dot denotes the derivative with respect to the affine parameter p. The 
second geodesic is defined by using the Christoffel connection of the given metric instead 
of the independent (metric) connection. This second geodesic minimizes locally the 'arc 
length'. Both geodesies coincide if the torsion only has a completely antisymmetric part, 
V = 0, M = 0. 



3 Vector fields preserving metric and connection 

Let f be a diffeomorphism with Jacobian matrix 

A%{x):= ~^{x). (18) 

If g^ v {x) is the metric tensor of a metric g with respect to the coordinates x^, then by 
definition (p is a (local) isometry if 

g, v {x) = (A T )/(x) g-^{x)) A\(x) . (19) 

Upon linearisation <p(x) = x + £(x) + o(£ 2 ), where £ = £ a d/dx a is a vector field, equation 



(19) becomes the Killing equation: 



Likewise if T x fiu (x) are the components of a connection (not necessarily the Christoffel 
connection of the metric) with respect to the coordinates x M , then by definition (p preserves 
the connection if 



T\ u (x) = (A- 1T )- X \x) A%(x) A\{x) T\^(x)) - (A T )J(x) — (A~ 1T ) p \x) . (21) 
Linearising we obtain: 



r _^_ T \ _ + ^ r A_ , g^pA i d 2 £ A =Q (22) 

* dx° ^ dx* ^ dx» » v + dx v ^ dx»dx» 1 } 

4 The most general Riemann-Cartan space invariant 
under 0(3) x M 3 

In order to combine Einstein-Cartan's theory with the cosmological principle we have 



to solve the Kiling equation (20), its analogue for the connection (22) and the metricity 
condition (16) for all vector fields generating the maximal isometry group of space. To 
simplify the calculations - and not because we particularly believe in flat spaces - we will 
take this group to be the 3-dimensional Euclidean group 0(3) k M. 3 generated by three 
infinitesimal rotations £ = yd/dz — zd/dy, zd/dx — xd/dz and xd/dy — yd/dx, and by 
three infinitesimal translations £ = d/dx, d/dy, and d/dz, where x, y, z are Cartesian 
coordinates. 
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The most general solution of the Killing equation is well known dr 2 = h dt 2 — a 2 (dx 2 + 
dy 2 + dz 2 ), with two positive functions a(t) and h(t). By a redefinition of the time 
coordinate t, we set h = 1. 



Equation (22) is the analogue of the Killing equation and describes connections in- 
variant under a vector field £. Technically it consists of 64 differential equations for each 
vector field, of which some are empty, some are identical, some are algebraic. Our task 



therefore is straight-forward and tedious: write down and solve the 6 x 64 equations (22) 
coming from the six vector fields £. For the three translations we get immediately that all 
64 components of V are independent of x, y and z. For the three rotations, it is sufficient 



to solve equation (22) for two infinitesimal rotations. It is then automatically satisfied for 



the third rotation because this rotation is the commutator of the two others. It remains 



to treat 2 x 64 equations. Finally, the most general connection solving equation (22) with 
the six vector fields has the following non-vanishing components: 

T t tt = A(t) } r* - = o(t)6(t)^, r ]t = c(t)5 iv (23) 
r tj = d{t)5, tj) r jk = f(t)e ljk , e 123 = l, (24) 



with five arbitrary functions of t: A, b, c, d, f. The metricity condition (16) yields with 
' :=d/dt, 

a' , b , „ 

A = Q, c = — d= -. (25) 
a a 

So far we have only taken into account transformations close to the identity and we still 
have to discuss space inversion. Since a and b are parity-even and / is parity odd, we 
must set / = 0. 

This connection reduces to the Christoffel connection if b = a'. 



We use the gauge transformation (15) with e a M = diagonal(l, a, a, a) to compute the 



non-vanishing holonomic components U) a b^ of the connection: 

= u\j = b Sij . (26) 
The Riemann tensor has the following non-vanishing components: 

. y i b 2 

R iOk = R l Q0k = — Sik, R l jki = ^ (^ikSjl — Sigdjk)- (27) 



The Einstein tensor has: 



b 2 ( b' b 2 



Goo = 3 - , Gij = - 2 - + - 5 ir (28) 
a z \ a a 1 J 

Note that it is symmetric but not covariantly conserved. 
The torsion tensor has: 

T\ = ^ . (29) 

In particular the vector part has only a time component, Vq = (b — a') /a, and the anti- 
symmetric and mixed parts vanish, A = 0, M = 0. This result agrees with the torsion 
found in references JTUJ 0] for spacetimes with maximally symmetric 3-spaces. 
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5 Hubble diagram 



To make contact with physics, we must solve the geodesic equations for co-moving galaxies 
and for photons [11]. For both, torsion decouples and they reduce to the geodesic equations 
with the Christoffel connection of the metric. Consequently the redshift is still given by 
z = clq/ a(t) — 1 with look-back time t and the apparent luminosity i is still related to the 
absolute luminosity of the standard candle L by 

'(*) = a 2 L m2 • (30) 

We have put the earth at the origin of the Cartesian coordinates and the supernova on 
the x-axis: 

x(t):= -jr. (31) 

It is remarkable, that the Hubble diagram in absence of peculiar velocities is entirely 
determined by the metric and does not feel torsion. This is not true if the standard 
candle has non- vanishing peculiar velocity and / or if the messengers are massive neutrinos 
instead of photons. 

On the other hand, the Einstein equations in presence of half-integer spin do feel 
torsion. Note however that the link between the Hubble constant Hq and d(z 2 £)/dz(0) 
is purely kinematical. It does not depend on the dynamics, i.e. the Einstein equations. 
Therefore it does not depend on torsion. This fact will be crucial to identify consistently 
the initial conditions of the Einstein equations. 



6 Matter, equations of state and field equations 

Before we can write down the Einstein equations and see how they are modified by torsion, 
we must talk about their right-hand sides. 

The most general energy- momentum tensor invariant under 0(3) x 1R 3 must sat- 



isfy an equation similar to the Killing equation (20). It therefore contains two arbi- 
trary functions of cosmic time, p(t) and p(t) interpreted as energy density and pressure: 
t^v = di&g(p,a 2 p,a 2 p,a 2 p). Note that we have already used the freedom of coordinate 
transformations to get rid of one function, h(t), in the metric. We are therefore stuck now 
with two functions in the energy-momentum tensor. But then we have three functions, 
a, p, p and only two independent components in the Einstein tensor and must get rid of 
one function. This is usually done by an equation of state pit) =: w pit). To simplify 
even more, one thinks of highly diluted particles and sets the pressure to zero, w = 0. 

Likewise the most general 0(3) x M 3 -invariant spin tensor has one function of time s(t) 
in the time component of the vector part, s(t) := ^SoabV ab - In Einstein-Cartan's cosmology 
for dust we then have four functions, a, b, p, s but only three equations, two of Einstein 
and one of Cartan. Therefore we need another equation of state. Our motivation is to do 
without dark matter. Visible matter consists mostly of protons, neutrons, electrons and 
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neutrinos, all of spin | and all with small densities today. Therefore we keep neglecting 
the pressure. We also neglect spin-spin interactions and assume the spin density to be 
proportional to the particle density or energy density: 

s(t) =: w s p(t). (32) 

Note that, while w is dimensionless, w s has units of time. For a single proton, we would 
expect w s ~ h/(mpc 2 ) ~ 10 -25 s. For several protons with some spins anti-parallel, w s 
could be even smaller. For a neutrino, w s should be larger. 

We have not tried to derive the spin density from any microscopic model. Kopczyhski 
[12] assumed a Weyssenhoff fluid which has the draw-back of producing a spin-density 
incompatible with maximal space-symmetry. Bohmer & Burnett [13] studied Elko spinors 
or "dark" spinors as microscopic model. These spinors are CPT-odd and couple to 
ordinary matter only via gravity and via the Higgs mechanism, motivating the second 
terminology. Furthermore, they produce a spin-density compatible with the cosmological 
principle. Poplawski [T3] uses a particular form for the spin density, which he derives from 
quarks and leptons. 

Now we are in position to write down the generalised Friedmann equations, i.e. the tt 
and the xx components of Einstein's equations, and Cartan's equation: 

= A + 8rrGp, (33) 
= A, (34) 
= 8nGw s p. (35) 







b 2 




3 








a~ 2 


V 




b 2 




+ 




a 




a 2 




a! - 


- b 


3 







We have three equations for three unknown functions, a, b and p. The first equation (33) 



is algebraic, the other two equations, (34) and (35), are first order differential equations. 
We use the algebraic one to eliminate p from Cartan's equation (35). Then we have 
a unique solution with two inital conditions a(0) = ao and 6(0) = 6q. We therefore 
have four parameters, a , b , A and w s . (We assume Newton's constant known.) These 



four parameters are not independent. With equation (33) we now eliminate b , Cartan's 



equation then relates czq, A and w s to the Hubble constant H . The remaining equation 



(34) allows us to compute the decceleration parameter go- We do not spell out this 
complicated expression for go preferring to obtain the Hubble constant from the Hubble 
diagram at z = (without using any dynamics) and to obtain the cosmological constant 
and the state parameter w s from a fit of the entire Hubble diagram. Let us write the 
relation between the Hubble constant and a , A and w s : 

1 = ^m0 + ^ao + 2fi s o — ^s0> (36) 
using familiar dimensionless quantities: 

8nGp A rr^rGp 
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In particular, we see that the scale factor today a has dropped out. This is well-known for 
cosmology with vanishing spatial curvature and remains true in presence of non-vanishing 
torsion. 



7 Data analysis 

To confront Einstein-Cartan's theory with experimental data, we have used the so-called 
Union 2 sample [15] containing 577 type la Supernovae up to a redshift of 1.4. Our 
results have been obtained using the full covariance matrix including correlations and 
systematic errors. The magnitude evolution of supernovae with redshift is written as 
M(z) = m s — 2.51og£(z) where m s is a normalisation parameter fitted to the data and 



z) the apparent luminosity (30). As in standard cosmology, the magnitude evolution 



with redshift is Hubble-constant free. 



The apparent luminosity is computed by solving the differential equations (34) and 



(35) where the algebraic equation (33) has been used to eliminate p. This system of 2 
coupled differential equations is solved numerically by using the Runge-Kutta algorithm 
[To*] with an adapted step in time such that the equivalent redshift step is much smaller 
than the experimental redshift error: Az < 10~ 5 . 

The best fitted cosmology is obtained using the MINUIT [17] package for the iterative 
X 2 minimisation with: 

X 2 = AM T V~ 1 AM (38) 

where AM is the vector of differences between measured and expected magnitudes and 
V the full covariance matrix. 

Since correlations between parameters are expected to be high, we marginalised over 
unwanted parameters and compute errors or contours by using the frequentist prescrip- 
tion [18] and solving x 2 — Xmin + s 2 (s = 1 for a la error on a single variable) where 
Xmin i s the minimum x 2 over all unwanted parameters. The Einstein-Cartan cosmology 
fit is performed with 3 free parameters (m s , Q m , Q s ) while is a derived parameter 



obtained from the Friedmann-like equation (36). All of the following results are given 
with marginalisation over the nuisance parameter m s . 

Table 1 presents the results for the fit of Einstein-Cartan's theory and for comparison 
of the pure Einstein theory. To get the same number of free parameters, the fit for the 
pure Einstein theory does not use the flatness constraint. The minimum x 2 f° r both fits 
are very close which implies that both theories are compatible with supernovae data at 
the same statistical level. The preferred value for Q s is 0.12 compatible with zero at a 
level better than 1 sigma. The effect of torsion is to lower the preferred value of Q m to 
0.09 which is statistically compatible either with the baryon density of 0.046 ±0.003 or the 
total matter density of 0.27 ± 0.03 published by the WMAP collaboration [19]. However, 
the spin density can contribute to the dark matter energy density to a certain amount. 
On the contrary the cosmological constant is only mildly modified by torsion. 

Figure [2](b) shows probability contours in the (f2 TO ,f2 s ) plane. As already mentioned, 
the degeneracy between both parameters is very high. As we can see, for negative spin 
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n qc;+u.24 

U.OJ_ 02 g 




0. 


530.0 


Einstein-Cart an 


u - uy -0.07 


o.83±8:±g 


12 +u ' 02 

U ' iz -0.22 


530.4 



Table 1: Fit results (la errors) for Einstein and Einstein-Cartan theories. No flatness 
constraint is imposed in Einstein's theory. 

density, the matter density increases while the cosmological constant density decreases 
(Figure [3]). This suggests that a negative spin density may act as dark energy density. 

We test this hypothesis by setting Q A equal to zero in the Einstein-Cartan equations 
and we redo the fit with only m s and £l s as free parameters. In Figure |2](a) the green 
(lower) curve shows the result of the best fit for Einstein-Cartan theory fixing Q A to zero, 
while the red (upper) curve shows the best fit with free cosmological constant energy 
density. In the first case, the minimum \ 2 is equal to 1188 for a number of degrees of 
freedom of 575 while in the second case the minimum x 2 is 530.4 for a number of degrees 
of freedom of 574. These results rule out the hypothesis of torsion to replace dark energy. 
Let us mention the analysis of reference [20], which contains a similar analysis using the 
parity-odd torsion [21]. They find that this part of torsion can mimic a cosmological 
constant. 

In Figure [3] we show the probability contours in the (fi m , ^a) plane for both Einstein- 
Cartan and Einstein theories. The degeneracy between both theories is almost orthogonal. 
The main reason is coming from the Friedmann-like equation for Einstein-Cartan theory. 
This constraint mimics an approximate flatness fl m + Q\ = 1. We notice that both 
contours cross at a confidence level of 39% at the usual ACDM point, corresponding to 
vanishing torsion. From a purely statistical point of view, the actual supernovae data do 
not allow us to discriminate between both theories. 

8 Conclusion 

It is well known that curvature is not an alternative to dark matter in the Hubble diagram, 
see Figure 3. On the other hand, within the error bars of the supernova luminosities, 
torsion a la Einstein-Cartan allows us to fit the data with a matter contribution of Q m = 
5 % corresponding to visible matter only. However the best fit has Q m = 9 % implying 
still a few percent of dark matter. More disturbingly, its corresponding state parameter is 
w s ~ Hq 1 ~ 10 17 s, 42 orders of magnitude away from the naive value w s ~ h/{m P c 2 ) ~ 
10- 25 s. 

To conclude, we think that the present result is encouraging enough to reconsider the 
rotation curves in galaxies and the CMB anlysis in the light of the Einstein-Cartan theory. 

Acknowledgement: It is a pleasure to thank Sami ZouZou for his advice. 
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Figure 2: (a) Fit results using the Union 2 Supernovae sample. The red (upper) curve corre- 
sponds to the Einstein-Cartan 3-fit (m s ,£l m ,£l s ) while the green (lower) curve represents the 
2-fit assuming a vanishing cosmological constant, (b) 39%, 68% and 95% confidence level contour 
in the (O m ,f2 s ) plane for the Einstein-Cartan 3-fit. 
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Figure 3: 39%, 68% and 95% confidence level contours in the (Q m ,fl\) plane for the Einstein- 
Cartan 3-fit (full contours) and for the Einstein 3- fit (full colored (elliptical) lines). The black 
(straight) line corresponds to the usual Einstein flatness. 
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